Classical and quantum chaos in a circular billiard with a straight cut 
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We study classical and quantum dynamics of a particle in a circular billiard with a straight 
cut. This system can be integrable, nonintegrable with soft chaos, or nonintegrable with 
hard chaos, as we vary the size of the cut. We use a quantum web to show differences in the 
quantum manifestations of classical chaos for these three different regimes. 
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I. Introduction 

In recent years, the dynamics of noninteracting par- 
ticles in two-dimensional (2D) billiards has been studied 
in many different billiard shapes [1-6]. In this paper, we 
study the classical and quantum dynamics of a particle 
in a closed circular billiard with a straight cut (we name 
this system the "Moon" billiard; see Fig. 1). The classi- 
cal system can exhibit both integrable and nonintegrable 
behavior. It also shows two distinct types of chaotic be- 
havior, both "hard chaos" and "soft chaos" as we change 
the size of the cut. Systems whose dynamics consists of 
a mixture of stable and unstable periodic orbits are said 
to exhibit soft chaos. If all periodic orbits in a system are 
unstable throughout, the behavior is called hard chaos [7]. 

The quantum version of these billiards has a discrete 
energy spectrum, and chaos (or nonintegrability) mani- 
fests itself in the quantum billiard in indirect ways. The 
statistics of energy levels for classically integrable systems 
is different from that for classically chaotic systems [2- 
5,8-11], and one commonly uses the spectral statistics of 
energy levels as signature of an underlying chaos. Also 
discrete symmetries of the system should be handled care- 
fully. 

In this paper, we focus on a diagnostic tool which 
has proven successful for spin systems. We calculate the 
quantum web [11-14] using about a hundred lowest-energy 
eigenstates, and observe patterns for classically different 
cases. The lattice-structure, which appears for classically 
integrable cases, breaks in different ways for cases with 
soft chaos and hard chaos. The Husimi function [15,16], 
which extracts the quantum Poincare section from a quan- 
tum state, is also calculated for some selected energy 
eigenstates to examine the quantum web more closely. 



with w. For example, Mi is a full circle and M\ is a half 
circle, and in general w is in the range of < w < 2. The 
classical Hamiltonian of a particle inside M w , using polar 
coordinates (r, 6), is 



H w (r,6,p r ,pe) 



Pr | Pi 

2m 2mr 2 



+ V*(r,6), (1) 



where V^(r,9) represents the hard- wall potential of the 
billiard with radius R and W = wR. To begin, we will 
study the full circle (M2) and the half circle (Mi). Both 
of cases are integrable. Then nonintegrable cases will be 
examined. 

A. Full-Circle Billiard, M 2 

The dynamics of a particle inside M2 is integrable be- 
cause there are two constants of motion, the energy E 
and the angular momentum pg. Given E and pg, the or- 
bit lies on a 2D torus (2-torus) in phase space. There ex- 
ists a canonical transformation to action-angle variables, 
{J r , Jg,4> r ,4>e), where actions are given by 



Je = 



R 2 



n 

2mE 



be I 
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Pe, 



(2) 

with the new Hamiltonian H' = E(J r ,Jg). We can also 
find the angle variables, but it should be noted that <f>g is 
not equal to 9. We can find explicit expressions of angular 
frequencies 4> r = io r and <pg — Lug as functions of E and 
pg using Eq. (2), 



OJ = — 



dE 

dJ r 



dJ r 
~dE 



2ttE 
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(3) 



II. The Classical Billiard 

In this section, we focus on the classical dynamics of 
the Moon billiard. We introduce a dimensionless param- 
eter w = W/R to characterize the system where W is 
the width of the billiard and R is the radius (see Fig. 1). 
Then we define M w as a circular disk with a straight cut 
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It is useful to look at periodic orbits (orbits for which 
the ratios of two angular frequencies are rational). We 
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classify periodic orbits using the notation (to, n) where m 
and n are relatively prime positive integers defined by 



I cog- 1 ( -Je^LJ] = ™ 



(5) 



where 2m < n [see Fig. 2(a)]. If we have a periodic orbit 
classified as (to, n), an infinite number of rotated periodic 
orbits also belongs to (to, n). Therefore periodic orbits in 
the full circle are non-isolated [1]. For any periodic orbit 
classified as (to, n), there are n bounces and to revolutions 
in one cycle. On the other hand, a non-periodic orbit will 
not come back to the starting point, and eventually fill 
the whole 2-torus. 

B. Half-Circle Billiard, M 1 

For the half circle, we still have two constants of mo- 
tion, E and p$. The range of 9 is reduced to a half, 
— -| < 9 < Tj, but we can still construct 2-tori on which 
orbits lie. For any orbit in Mi there is a corresponding 
orbit in M 2 . (If there is an orbit in M 2 , folding M 2 in half 
gives us a corresponding orbit in Mi.) Hence we can use 
the results found for M 2 to describe some periodic orbits. 
A periodic orbit is classified as (to, n)' when the corre- 
sponding orbit in M 2 is (to, n). Unlike for M 2 , periodic 
orbits in M\ are isolated, since there is no rotational sym- 
metry. In the group (to, n)' , there are an infinite number 
of different periodic orbits, but only a few play an impor- 
tant role when w is slightly less than 1. In Fig. 2(b), those 
periodic orbits are shown. The stabilities of these peri- 
odic orbits are all neutral (neither stable nor unstable) 
like those in M 2 . 

It is useful to see periodic orbits in the space of E 
and \pe\. A point in this space corresponds to a group of 
orbits with constants of motion E and pg. When Eq. (5) 
is satisfied, the group consists of periodic orbits. Thus 
the condition to have periodic orbits is given by Eq. (5), 
and is plotted in Fig. 3. They are densely populated in 
the classically allowed region. The classically forbidden 
region in this space is given by E — pg/(2mR 2 ) < 0. In 
integrable cases, orbits are a point (zero-dimensional) in 
(E, \pg |)-space since E and p@ are conserved, but as the ro- 
tational symmetry breaks by changing w, p 2 e is no longer 
conserved and orbits become one-dimensional. We will 
see later that Eq. (5) gives us the condition for primary 
resonances when w is slightly less than 1. 

C. M w when < ro< 1 or 1 < w< 2 

Except two cases studied above (w = 1, 2), the system 
is nonintegrable, because E is an only constant of motion 
in the system of two degrees of freedom. The stabilities of 
periodic orbits in M w give us one way to understand the 
dynamical behavior of the system. The simplest periodic 
orbit for any w is an orbit with two bounces, going back 
and forth [for example, (1,2)' in Fig. 2(b)]. To calculate 



the stability of this orbit, we need to imagine a new bil- 
liard M 2 , which is a composite of circular parts of two 
M w 's facing each other. It is easy to see that Mf is just 
Mi. Orbits in M w have the correspondence with orbits 
in M 2 as orbits in Mi correspond to orbits in M 2 = Mf . 
The stability of this two-bounce orbit in M 2 can be cal- 



culated from a 2 x 2-matrix m, acting on (66, 
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on the boundary, which decides the deviation after two 
bounces [1] 



2(2w- l) 2 
4(2io- 1)(1 - 
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Aw(l - 2w) 
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(6) 



The eigenvalues of m, A± , are given in terms of the trace 
of m, 

A± = i{Trm± [(Trm) 2 -4] 1/2 }, (7) 

where we used det|m| = 1 since m is area-preserving. 
The orbit is neutral when |Trm| = 2, stable when 
|Trm| < 2, and unstable when |Trm| > 2. The two- 
bounce orbit is neutral when w = 0.5, w = 1, or w = 2, 
stable when < w < 0.5 or 0.5 < w < 1, and unstable 
when 1 < w < 2. We have seen that all periodic or- 
bits are neutral when the system is integrable (w = 1,2) 
but w = 0.5 is a special case as we will see later. From 
this result, we can predict that there are no stable pe- 
riodic orbits in the billiard Mi <tu<2 , and periodic orbits 
in M <to<i are either stable or unstable, except when 
w = 0.5. (Ergodicity of the billiard M\ <w< i has been 
proven by Bunimovich [17].) The system shows hard 
chaos when 1 < w < 2, and soft chaos when < w < 1. 

The Poincare surface of section (PSS) is one way 
to observe the chaos. Here PSS is a two-dimensional 
surface (9,pg) at r = R along the circular boundary 
{-6 ma x <9 < 9 max , -RV2mE < p e < R^j2mE). Each 
orbit gives a point in (9 , pg)-sp&ce every time it touches 
this surface. Therefore PSS becomes a 2D area-preserving 
map. In fully chaotic (ergodic) cases, points generated 
by an orbit will fill out the whole allowed space in the 
PSS. In cases of soft chaos, we observe some structure. 
Some orbits generate island chains and some orbits will 
fill some remaining regions in a chaotic manner. Since the 
KAM theorem [7,11] does not apply in this system due to 
singularities of the boundary, the existence of KAM tori 
seperating island chains is not guaranteed even for small 
S when w = 1 — 5. 

In Fig. 4, we show PSS's for various w's. In Fig. 4(a), 
one orbit is filling the whole region when w — 1.01. This 
is a fully chaotic case. Figure 4(b) is an integrable case 
when w = 1. When w = 0.99 [Fig. 4(c)], we see is- 
land chains centered at the positions of periodic orbits 
along with stochastic diffusion in the remaining region. 
The resonance condition, Eq. (5), gives us locations of 
island chains. For example, (1,2)' gives us the loca- 
tion p e /(RV2mE) = 0, (1,4)' gives us p e /{Ryj2mE) = 
±1/Vt, (1,3)' g ives u s p e /{RV2m~E) = ±0.5, and (1,6)' 
gives us pe/(R^j2mE) = ±\/3/2, and so on. As w de- 
creases [Figs. 4(d) and 4(e)], we see increasing chaotic 



2 



region and some remaining island chains. When w — 0.5 
[Fig. 4(f)], periodic orbits that have neutral stabilities re- 
side in two axes, (9 = 0)-axis and {pg = 0)-axis. One 
orbit starting near the center is filling out almost all re- 
gion. This is the most chaotic case when < w < 1. As 
w decreases further, we see bigger regular regions in PSS 
in Figs. 4(g) and 4(h). 

III. The Quantum Billiard 

In this section, we study the quantized M^-billiard us- 
ing a quantum- web analysis with about a hundred lowest- 
energy eigenstates, and we also look at some individual 
energy eigenstates using Husimi plots. 

A. Numerical Method 

The Schrodinger equation for this 2D M^-billiard is 
the Hclmholtz equation, 

(V 2 + fc 2 )*(r) = 0, (8) 

with the Dirichlct boundary condition, ^ = 0, on the 
boundary of M w , B w = dM w , where k 2 — 2mE/h 2 and 
V 2 = ^2+i^ + ^^2 using polar coordinates. 

For the classically integrable cases, this equation can 
be solved analytically. The Hamiltonian H w (w = 1 or 
2) and the angular momentum pg (p$ for the half circle) 
commute. They are generators of continuous symmetry 
transformations, the time translation and the rotation. 
For a full circle Mi, we can find energy eigenstates which 
are simultaneous eigenstates of H2 and pg, 



(r|i,fc)(/)«J l (^) 



(9) 



where k is a positive integer, I is an integer, and an- is 
the fcth zero of the the Bessel function Ji(x). And en- 
ergy levels are given by E^P — h 2 af k /(2mR 2 ). There 
exist two-fold degeneracies when I ^ since the system 
also has the parity symmetry and [pg,P\ ^ (P is the 
parity operator). Then wc can find another set of energy 
eigenstates, simultaneous eigenstates of H 2 , pe 2 ' , and P, 



(r\l,k,+)W oc J t (<*nf) C0S (W) (Z>0), 
(r|Z,fc,-)(/) cx J, ("if ) sin(Z0) (I > 1), 



(10) 



where — tt < 9 < it. The latter will be used in the calcu- 
lation of the quantum web. 

For a half circle Mi, energy levels are the same as 
those of the full circle without levels with I = 0, and the 
energy eigenfunctions are 



(r|J,fc,+)W oc J x (suj:) cos(W) (/ = 1,3,5,...), 
(r|i,fe,-)W cx Ji (^jr) sin(Z#) (Z = 2,4,6,...), 



Eq. (8) must be solved numerically. We use the boundary 
element method (BEM) [2,18-20] to solve this nonseper- 
able 2D partial differential equation. It is an efficient 
way to solve boundary-value problems, because in BEM 
a 2D equation with boundary condition becomes an in- 
tegral equation in one variable along the boundary. The 
method we use is briefly outlined below. 

We use a Green's function, G(r,r') — — - wWf 
r'|), which satisfies 



(11) 



where — ir/2 < 9 < ir/2. Here there is no degeneracy. 

For classically nonintegrable cases, pg 2 no longer com- 
mutes with H w (but still [P,H W ] = for any w). Here 
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(V 2 + fc 2 )G(r,r') = 5(r-r'), 



(12) 



where H> [x) is the Hankel function of the first kind. 
We multiply Eq. (8) by G(r, r'), and multiply Eq. (12) by 
\&(r'). After subtracting one from the other, integrating 
over the area of M w with respect to r', and using Green's 
theorem, we finally get 



ds' G{r, r') u(s') = 



*(r), 

£*(r), 
0, 



r inside B w , 

r on B w , 

r outside B w , 



(13) 



where s' is arc-length along B w = dM w , u(s') = n' • 
V'4'(r(s')), and n' is the outward normal unit vector to 
B w at r'. With r lying on B w , by taking the normal 
derivative n • V on both sides of Eq. (13), we obtain 



u(s) 



ds'u(s') (n- V)G(r, r') 



(14) 



One way to solve this equation is discretizing B w by di- 
viding it into N segments. Then Eq. (14) becomes the 
matrix equation A • x = x, where A = A(fc) is an 
(N x A r )-matrix and x is an A-component vector rep- 
resenting {u(sj)|l < i < N}. For given w, energy levels 
of the system, E n = h 2 k 2 /2m (n > 1), can be found by 
solving the equation, det |A(fc) — I| = 0. For each energy 
level E n , we can obtain {u n {si)} by finding an eigenvec- 
tor of A(k n ) with a near-zero eigenvalue. Since the nu- 
merically obtained energy levels in this way always have 
some uncertainty, "degeneracy" (which is actually near- 
degeneracy) can occur when the difference between two 
adjacent exact energy levels is less than the uncertainty. 
In these cases, we can find two sets of {u(si)} with two 
near-zero eigenvalues of A. Therefore, looking at eigen- 
values of A(k n ) can be an easy way to check numerically 
for near-degeneracies of an energy level, E n . 

For given w, we found energy levels {E n \n > 1} and 
normal derivatives, on the boundary, of corresponding en- 
ergy eigenfunctions {^> n (r) = (r\n)} where |n)'s are en- 
ergy eigenstates. Then from Eq. (13) we can calculate the 
energy eigenfunction inside, 

* n (r) = j/ ds'H^(k n \r-r'\)u n (s'). (15) 
4 Jb. 
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Using Eq. (15), we can also calculate (r\p"g \n), which will 
be used in calculations of quantum webs in the next sub- 
section, 



(r\P9 2 \n) 



ft 2 J^n(r) 
ds 



4 IB 



^H^(k n \r-r'\)\ u n (s'), 
(16) 



H^\k\v-v'\) 



where 

+k 2 [^|r-r'|] 2 ^ 1) "(fc|r-r'|) 

(17) 

and we use 



|r — r'| = \J (r cos 9 — x') 2 + (r sin 9 — y 



A2 



(18) 



and 



H^\x) 



(19) 



In calculations of and -^^ n , care must be taken when 
r is close to the boundary because H^\\r — r'|) diverges 
as |r — r'| goes to zero. 

B. The Quantum Web 

The quantum-web analysis can be used to observe the 
manifestations of classical chaos in quantum mechanics, 
and until now has been used primarily in spin systems 
[13,14]. Here we will look at three different cases: clas- 
sically integrable cases (w = 1,2), nonintegrable cases 
showing soft chaos (0 < w < 1), and nonintegrable cases 
with hard chaos (1 < w < 2). 

For classically integrable cases (w = 1,2), we have 
seen in Sec. II that there are two constants of motion, 
E and \pg\, and that we can find two action variables 
(J r ,Jg). (There exists a nonlinear map from (J r ,Jg)- 
spacc to (E, |pe|)-space.) In Fig. 3, a classical orbit ap- 
pears as a point in (E, |pg|)-space. In quantum mechanics, 
there exist simultaneous eigenstates of two operators H w 
(w = 1 or 2) andpg 2 [see Eqs. (10) and (11)]. We can con- 
struct a 2D space, in which a pair of eigenvalues (Eik, \lfi\) 
of each eigenstate is plotted point. In Figs. 5(e) 
and 5(1), we observe the structure of a deformed lattice. 
This can be understood as a mapping from an almost per- 
fect 2D lattice in ( J r , Jg)-space to the deformed lattice in 
(E, |pe|)-space. This almost perfect lattice-structure can 
be explained from Einstcin-Brillouin-Kcllcr (EBK) semi- 
classical quantization, 

J r ~ (n r — (n r : positive integer), 

J e = p g = ITi (I : integer), (20) 

where J r is a very good approximation and Jg is exact. 
This is the quantum analogue of Fig. 3. Each point in 



Fig. 3 represents a 2-torus. Thus, we can find a 2-torus in 
classical phase space corresponding to a quantum eigen- 
state, and then each eigenstate here corresponds to a 
group of orbits that are on this 2-torus. 

For classically nonintegrable cases, [pg 2 , H w ] ^ when 
w 7^ 1 or 2. However, we can still calculate an expecta- 
tion value of pg 2 for an energy eigenstate |n), (n\p g 2 \n), 
numerically using Eqs. (15) and (16), 



(n\pg 2 \n) 



d 2 r (n\v){v\pg 2 \n) 



s 7, (2i) 



These values can be interpreted as time-averages of pg 
[12]. When there is an accidental degeneracy (or near- 
degeneracy), we find expectation values from eigenval- 
ues of the matrix representation of pg in the degener- 
ate subspace. In this way, we obtain a pair of values 



(E n , y (n\pg \n}) for each energy eigenstate. These points 
can also be plotted in a 2D space as a quantum web. We 
expect that the lattice-structure that exists for integrable 
cases will be broken because EBK quantization doesn't 
apply to nonintegrable cases. 

In Figs. 5(a)-5(d), quantum webs are shown for cases 
of soft chaos. When w — 0.99 [Fig. 5(d)], we observe 
breaking of the web near conditions of primary resonances 
in classical mechanics. We see patterns of crossing near 
(1,4)'- and (1, 3)'-resonances. Although "regular" parts 
still exist, some layers seem to start to shift near resonance 
conditions. We can roughly estimate the energy value at 
which the effect starts to be seen for each resonance con- 
dition by measuring the width of island chains, Apg, in 
Fig. 4(d). For example, (1, 2) '-resonance has the biggest 
width, and next (1, 4)'-resonance, and so on. Because pg 
is scaled by (RV^mE)^ 1 in Fig. 4, Apg is proportional 
to \[E. When Apg(— AJg) > K, the resonance can be 
clearly seen in the quantum system, and we can roughly 
obtain an estimate of the minimum energy at which each 
resonance is in effect. The smaller the width of an island 
chain and the lower the energy, the less likely to find the 
web broken near the curve of the particular resonance. 
When w = 0.7 [Fig. 5(c)], we see similar patterns as in 
Fig. 5(d). When w = 0.5 [Fig. 5(b)], the classical system 
has a large chaotic region in the PSS, and has periodic 
orbits with neutral stability, which reside on two axes, 
6 = and pg = 0. The quantum web, however, is quite 
regular although the structure looks different from those 
of integrable cases. It looks more like a structure of lay- 
ers. When w = 0.3 [Fig. 5(a)], the web is similar to that 
of Fig. 5(b). 

In Figs. 5(f)-(k), we show quantum webs for cases 
of hard chaos. When w = 1.01 [Fig. 5(f)], the lattice- 
structure is still intact except a little kink, although this 
is the fully chaotic case classically [see Fig. 4(a)]. We 
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observe the lattice-structure quickly collapses as we in- 
crease w. When w = 1.5 [Fig. 5(i)], the structure is very 
irregular except for four regularly placed points near the 
top-right corner. (Some of eigenstates noted by arrows 
here will be examined in the next subsection using Husimi 
plots.) The case of w = 1.7 [Fig. 5(j)] is the most irreg- 
ular quantum web among cases shown. When w — 1.9 
[Fig. 5(k)], we observe splitting of degeneracies and also 
quick collapse of lattice-structure from an integrable case 
uu = 2 [Fig. 5(1)]. As we have seen so far, the lattice 
structure tends to collapse quickly irrespective of energy 
in cases of hard chaos, but there also exist small remnants 
of regularity in some cases. 

C. Quantum Poincare Section 

The Husimi plot provides a method of extracting a 
quantum Poincare surface of section (QPS) from a quan- 
tum state [15]. The QPS is a quantum analogue of PSS, 
which we have obtained in Sec. II. The Husimi function 
of an ID system corresponding to a state |\&) is defined 
as 



F(x ,Po) = \(x ,po\*)\' 



(22) 



where |xo,po) is a coherent state with a representation in 
configuration space, 



(x\x ,p ) = 2 



1 



7T(T Z 



exp 



(x - XqY 



+ jVo{x - x ) 



(23) 

In 2D billiards, Eq. (22) can be modified to create a 
Husimi function using the coordinate along the bound- 
ary [16]. For example, for M^-billiard along the circular 
part of the boundary, the Husimi function is defined as 



F n (6 ,peo) 



<W (9 , P90 \9')u n (9') 



(24) 



wh ere - 6 max < 6 < max and -Rs/2mE < p eo < 
RV2mE. Here (9'\9o,peo) has the same form as in 
Eq. (23) with a given by the value [9 max fi/ 1 (Rs/2mE)] 1 / 2 , 
and u n (9) is the normal derivative of the energy eigenfunc- 
tion on the circular part of the boundary, J^„(r)| r= ^. 

In Fig. 6, we show Husimi plots for selected eigenstates 
for three cases (w = 0.5, 0.9, 1.5). And, in Fig. 7, we show 
probability densities of wave functions, |^I'„(r)| 2 , for some 
of eigenstates chosen from Fig. 6. (Most eigenstates cho- 
sen in Fig. 6 can be found in Figs. 5(b), 5(d), and 5(i) with 
arrows pointed to them.) 

Figures 6(a)-6(d) show Husimi plots of energy eigen- 
states when w = 0.5. The pattern of PSS shown in 
Fig. 4(f) can be seen in these plots. We can get some 
information on chosen eigenstates from the quantum web 
[Fig. 5(b)]. The eigenstate for Fig. 6(a) is on the outer 
part, and the eigenstate for Fig. 6(b) is on the inner part 
of the quantum web. Figure 6(a) shows chaotic region of 
PSS, and Fig. 6(b) seems to correspond with a two-bounce 



orbit with neutral stability, which we can observe in the 
wave function [Fig. 7(a)]. The eigenstate for Fig. 6(c) is in 
the middle of the quantum web, and its Husimi plot and 
wave function [Fig. 7(b)] lie between two extreme cases 
above. The eigenstate for Fig. 6(d) has a relatively high 
energy, but the structure is similar to Fig. 6(c). 

Figures 6(e)-6(h) show Husimi plots of energy eigen- 
states when w — 0.9, and each eigenstate is picking 
up a classical primary resonance shown as an island 
chain in PSS [Fig. 4(d)]. The eigenstate for Fig. 6(e), 
located near the (1, 2) '-resonance in the quantum web 
[Fig. 5(d)], shows the pattern of the island chain of the 
(1, 2) '-resonance, although it is only the 9th highest en- 
ergy eigenstate. The eigenstates of Figs. 6(f) and 6(g), 
located at the crossing of two layers near the (1,4)'- 
resonance, show the pattern of the island chain of the 
(1, 4)'-resonance. Wave functions [Figs. 7(c) and 7(d)] of 
these states show the trace of unstable and stable peri- 
odic orbits [see Fig. 2(b)], respectively. The eigenstate 
for Fig. 6(h), which is on the (1, 3)'-resonance, also shows 
the pattern of the island chain of the (1, 3)'-resonance. As 
expected earlier, the (1, 2)'-resonance is observed in the 
Husimi plot at lower-energy eigenstates than the (1,4)'- 
resonance. 

Figures 6(i)-6(l) show Husimi plots of energy eigen- 
states when w = 1.5. The eigenstate for Fig. 6(i), located 
in the inner part of the quantum web [Fig. 5(i)], seems 
to be picking up the periodic orbit with neutral stability, 
which can be seen clearly in the wave function [Fig. 7(e)]. 
The eigenstate of Fig. 6(k), which is one of four regularly 
placed points in the quantum web, shows a whispering 
gallery state [Fig. 7(f)]. All four of these regularly placed 
eigenstates show similar Husimi plots. A relatively high 
energy eigenstate for Fig. 6(1) shows a more uniformly dis- 
tributed Husimi plot like the corresponding PSS in clas- 
sical mechanics. 

IV. Conclusions 

The quantum web is the quantum representation of 
(E, |p#|)-space in Fig. 3. In the regime of soft chaos, we 
observe that the lattice-structure obtained for the inte- 
grable case starts to break near the primary resonance 
conditions obtained from classical mechanics as the width 
parameter, w, decreases from w — 1. The effect of reso- 
nances is greater when energy is higher and the width of 
an island chain in the PSS is greater. Even for the most 
chaotic case in the soft-chaos regime, layer-structure re- 
mains. In the hard-chaos regime, the regular quantum 
web collapses more quickly regardless of energy as the 
width parameter, w, increases from w = 1. 
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when R<W<2R 



when 0<W<R 



Figure 1: Geometry of the "Moon" billiard. When W = 
2R, it is a full circle, and when W = R, it is a half 
circle. Here 9 m ax is given by the equation, cos 9 m ax 

11 II': It. 



(a) 




(1,2) (1,3) (1,4) (1,5) (2,5) 
(b) 




(1,2)' (1,4)' (1,3)' (1,6)' 



Figure 2: Closed orbits in integrable cases, (a) In a full 
circle, periodic orbits can be classified as (m, n) where 
orbits have m bounces and n revolutions in a cycle, (b) 
In a half circle, we can use the notation of a full circle to 
classify periodic orbits as (m, n)'. 




Figure 3: For a half circle, the condition in (E,\pg\)- 
space to have periodic orbits (m,n)' is given by Eq. (5), 
and some are plotted. These curves are densely populated 
in classically allowed region. 
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Figure 4: The Poincare surface of section of M^-billiard 
varying w where w = W/R. We observe ergodic motions 
when 1 < w < 2, and generic chaotic behaviors when 
< w < 1. (a) Ergodic when w = 1.01 with one orbit, 
(b) Integrable when w = 1. (c) w — 0.99. (d) w = 0.9. (c) 
w = 0.7. (f) w = 0.5. Closed orbits residing in (9 = 0)- 
axis and (pe = 0)-axis have neutral stabilities. One orbit 
is filling almost all space except two axes, (g) w = 0.3. 
(h) w = 0.1. 



Figure 6: Husimi plots for given cigenstates. (e = 

E /(Sp)) ( a ) w = °- 5 and 6 = !43.28. (b) w = 0.5 
and e = 740.79. (c) w = 0.5 and e = 1156.07. (d) w = 0.5 
and e = 9992.23. (e) w = 0.9 and e = 62.513. (f) w = 0.9 
and e = 365.64. (g) w = 0.9 and e = 371.99. (h) w = 0.9 
and e = 617.98. (i) w = 1.5 and e = 258.03. (j) w = 1.5 
and e = 284.94. (k) w = 1.5 and e = 268.80. (1) w = 1.5 
and e = 842.69. [A square at the top right corner of each 
plot represents the size of h (Planck constant). All cigen- 
states except (d) and (1) are pointed in quantum webs in 
Fig. 5 by arrows.] 



Figure 5: The quantum web of M„,-billiard varying w. 
Even and odd parity eigenstates are shown together, (a) 
w = 0.3. (b) w = 0.5. (c) w = 0.7. (d) w = 0.9. The 
condition of primary resonances in classical mechanics are 
also shown, (e) w = 1.0. Classically integrable case, (f) 
w = 1.01. (g) w = 1.1. (h) w = 1.3. (i) w = 1.5. 
(j) w = 1.7. (k) w = 1.9. (1) w = 2.0. Classically 
integrable case. There are two-fold degeneracies when 
pe ^ 0. [Arrows in (b), (d), and (i) indicate states which 
will be studied using Husimi plots in Fig. 6.] 



Figure 7: The probability density of some energy eigen- 
states. (e = E/(^)) (a) w = 0.5 and e = 740.79. 
[See Fig. 6(b).] (b) w = 0.5 and e = 1156.07. [See 
Fig. 6(c).] (c) w = 0.9 and e = 365.64. [See Fig. 6(f).] (d) 
w = 0.9 and e = 371.99. [See Fig. 6(g).] (e) w = 1.5 and 
e = 258.03. [See Fig. 6(i).] (f) w = 1.5 and e = 268.80. 
[See Fig. 6(k).] 
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